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Introduction 

Consider a (scalar) ODE £ of order n > 1 in J n M (the Jet bundle of order n over 
the manifold M of independent and dependent variables); it is well known that if this 
admits a (Lie-point) symmetry, then it can be reduced to an equation of order (n — 1). 
More generally, if it admits a p-dimensional Lie group of symmetries G, then it can 
be reduced to an equation of order (n — q), with q < p depending on the algebraic 
structure of the operation of G and its Lie algebra Q [1, 6, 13, 23, 24, 27]. 

The approach to symmetry reduction for ODEs goes roughly speaking as follows: 
if £ admits the vector field X (on M) as a symmetry (generator), then £ is invariant 
under Y = X {n \ the prolongation of X to J n M. Thus £ can be expressed in terms 
of the differential invariants of Y. Or, this can be recursively generated starting from 
those of order zero and one thanks to the "invariants by differentiation" property (see 
[1, 6, 13, 23, 24, 27]). 

Using symmetry adapted coordinates, say w and y with y the independent 
variable, the equation will not depend explicitly on the coordinate along the vector 
field X, say w, and thus its order can be reduced by one passing to consider it as an 
equation for z = w y and its derivatives. 

It was observed by Muriel and Romero back in 2001 that the same scheme works 
if Y is not the (standard) prolongation of X but instead some kind of "twisted" 
prolongation, depending on an arbitrary smooth (by this we will always mean C°°) 
function A : J X M — > R and hence called "A-prolongation" [15, 16, 17]. 
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The key to this result is that A-prolonged vector fields still have the property of 
"invariants by differentiation" , so the scheme working in the case of invariance under 
standardly prolonged vector fields can still be applied. 

The same holds for systems of ODEs, except that here a symmetry will lead to 
reduction by one of the order of one of the equations; in the case of systems of n first 
order equation, a symmetry will lead to reduction to a system of n — 1 equations. For 
multi-generators symmetry group, as usual the attainable reduction will depend on 
the algebraic structure of the operation of the group. 

It was shown by Pucci and Saccomandi [25] that for scalar ODEs the A-prolonged 
vector fields are essentially the only ones sharing the "invariants by differentiation" 
property; they also provided a neat geometrical understanding of this fact. 

In the present note we want to generalize the Muriel-Romero approach, and 
consider cases where the equation (or system) is invariant under a set of vector fields 
Y(j) defined on J n M , and obtained from vector fields in M under a further modified 
version of A-prolongations. More precisely, in this case the modified prolongation 
operation will not act on the single vector field, but rather on the set of vector fields 
(that is, prolongation of each of them will involve the other ones). This "joint- A" 
prolongation will depend on a matrix a defined by a set of smooth functions <Jij 
on J X M, and will be therefore also denoted as a-prolongation (in the same way as 
A-prolongations took their name from the function A). 

We will show that even in this case the standard approach to reduction sketched 
above is still valid, and hence the equation can be reduced. Afore precisely, we will 
show that the "invariants by differentiation" property still holds. This will follow 
by explicit algebraic computation, but later on we will also discuss the geometrical 
meaning of our result, generalizing the result by Pucci and Saccomandi [25]. 

Finally, we note that the same construction with "joint A-prolongation" is 
considered in a companion paper [9]; however there the assumptions, in particular 
on the vector fields and the functions, arc different: the two papers are related, but 
deal with different situations and lead to different results. 

Acknowledgements. This work was started in the course of a visit of GG and SW 
in Pisa; they would like to thank the Universita di Pisa, and especially GC, for their 
hospitality in that occasion. 

1. Preliminaries and notation 

We will firstly recall some basic notion, also in order to set some general notation to 
be used in the following. 

1.1. Equations, solutions, symmetries 

We will only consider ordinary differential cquation(s); the independent variable will 
be denoted as x € R, the dependent one(s) as u € U = R or u a e U C R p in the multi- 
dimensional case. We denote by M — X x U the phase bundle, and by (J k M, Tr k , M), 
or J k M for short, the associated jet bundle of order k. 

A differential equation £ of order n is a map F : J n M — > R, and is naturally 
identified with the solution manifold S(£) = F' 1 ^) C J n M. In the case of l- 
dimcnsional systems £, we have I maps F^ : J n M — > R (or cquivalently a map F : 
J n M -> R e ) and a solution manifold S(£) = F^O) = (F 1 )- 1 ^) n ... n (-F^ _1 (0) C 
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J k M. Note that for p-dimensional dependent variables, J n M has dimension (np + 1), 
and a system of t independent equations identifies therefore a solution manifold of 
dimension (np + 1 — £ ) . 

A function / : X — > £7 is a solution to the differential equation(s) 5 under study 
if and only if its n-th prolongation lies entirely in the solution manifold S(£). 

Let us now consider a vector field Y on J n M; we say that £ is invariant under 
Y if and only if its solution manifold is; that is, Y : S(£) — > TS(£). This can also be 
cast as the condition [Y(£)} s ^ = 0. 

If Y is the prolongation of a (Lie-point) vector field X on M, Y = X^ n \ we say 
that X is a symmetry for £ (more precisely, this would be a symmetry generator; we 
will adopt this standard abuse of notation for ease of language). The condition for X 
to be a symmetry is therefore [X( n \£)] s , £ , = 0. 



1.2. Local coordinates 

We will consider local coordinates (x,u a ) in M = X x U, and correspondingly local 
coordinates {x,u%) (with k = 0, ...,n), where u% := (d k u a /dx k ), in J n M. 

A general vector field on J n M will be written in local coordinates (here and below 
we will use the Einstein summation convention; the notation denotes u together 
with its derivatives of order up to k) as 

y = cm£ + r k {xM k) )^-, a) 

this is the prolongation of 

X = ^,u)l + ^ X ,u)^ (2) 
if and only if the coefficients ipl satisfy the (standard) prolongation formula 

n+x = D *r k - <+i d x z ■ r a = ^ a . (3) 

The notation D x identifies the total derivative with respect to x, 

D x = d x + ut +1 A . (4 ) 



1.3. X-prolongations 

After the work of Muriel and Romero [15, 16] it is frequent to also consider A- 
symmetries of ODEs (see also [17, 18, 19, 20, 21] for their later work, and the 
bibliography given in [11]). 

A vector field Y written in local coordinates in the form (1) is a A-prolonged 
vector field if its coefficients satisfy 

r k+1 = o x r k - D *z + a (r k - (5) 

with A a smooth function A : J 1 M — > 77. We also say that Y is the n-th A-prolongation 
of X (see (2)) if ipfi = ip a i.e. if X is the restriction of Y to M. 

If the A-prolongation Y of X leaves an equation (or system) £ invariant, we say 
that A is a A-symmetry for £ . 
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1.4. The "invariants by differentiation" property 

When considering a vector field in J n M, it is of interest to know its invariants, i.e. 
the functions £ : J n M — > R such that Y(() = 0; these are also called differential 
invariants (to distinguish them from "geometrical invariants", i.e. functions defined 
on M rather than on J n M). 

For prolonged vector fields, it is well known that once we know differential 
invariants of order zero and one, say r] and (0 respectively, we can generate recursively 
differential invariants of any order. (Note that starting from differential invariants 
of order one will also be needed in the case of A-prolongations in order to take into 
account the properties of A.) In fact, the functions 

- <« 

turn out to be automatically invariant under Y if this is a prolongation and if 77 and 
(k were [1, 6, 13, 23, 24, 27]. This is also known as the invariants by differentiation 
property (IBDP). 

This property remains true in the case of vector fields which are not standard 
prolongations but are A-prolongations, and this fact is at the basis of the approach by 
Muriel and Romero; see e.g. their papers [15, 16] for an algebraic proof. 

The validity of IBDP was understood in geometrical terms by Pucci and 
Saccomandi [25]; their argument can be recast in the light of the approach by 
some of us [8] as follows. Any A-prolonged vector field is collinear to a standardly 
prolonged vector field (the standard and the A-prolongation being applied to different 
vector fields; the factor involved in the transformation can be nonlocal); however the 
invariants, and hence also the IBDP, only detect the direction of vector fields, not their 
magnitude (provided this is nonzero, of course), and hence are the same for collinear 
vector fields. 

In the following we will essentially generalize this remark to the case where we 
have several vector fields, so that the scope for variations in the generators of the 
resulting distribution is ampler. 

2. Joint prolongations 

We consider a set X = {X^, i = 1, r} of vector fields X^ in involution on M, i.e. 
such that 

[X (i) ,X {j) ] = 4(x)X k (7) 

with /x*. = —fiji smooth functions on M. We will also consider vector fields Y^ 
on J k M, which will be some kind (to be specified in a moment) of generalized 
prolongation of the X^ . 

In local coordinates, and with := (d k u a /dx k ), these will be written as 

d d 
X {i) - (x, u) — + ip a {i) (x, u) — , (8) 

% =£wM£ + ^i),fcfr« (fc) )fl^. ( 9 ) 

where we set ^ i) = ip a (i) . 

In the following we will also use the shorthand notation 

d k a := (d/dut) ; (10) 
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with this eq. (9) reads 

% - Z(i)dx + ^h k d k a , (11) 

Definition 1. Let a = {a\ = a^, i,j = 1, r} be a set of r 2 smooth real functions 
on J X M. The vector fields y = {Y^, i = 1, ...,r} on J k M (i = 1, ...,r), written in 
local coordinates as in (9), are said to be jointly X-prolonged (or a -prolonged) if the 
coefficients tp^ k satisfy, for k > 0, 

4>hk+i = {p**h,k - + °a {r U) , k - <+itu)) • (i2) 

If tpf^Q — ip%\, we say that the y = {Y^} are the joint X-prolongation (or a- 
prolongation for short) of the X = {X^}. 

In the following we use the notation with lower indices (cry) for ease of writing 
when no confusion arises, resorting to the one with upper and lower ones (a^) when 
it becomes convenient to keep fully track of covariant and contravariant indices. 

Remark 2.1. Note that the notion of cr-prolongation refers to a set of vector fields, 
not to a single one. 

Remark 2.2. As cry : J X M — »■ R, this notion represents indeed a generalization of 
A-symmetries. At the same time, the fact we introduce matrices rather than scalar 
functions to describe the "twisting" of the prolongation operation makes this similar 
to /i-symmetries (in their ODEs version) [11]. Remark 1.1 above guarantees we are 
considering a really new notion, and we will see below it gives new results. 

Lemma 1. // the y are a -prolonged, they satisfy 

[Y(i),D x ] = CTijVy) - (D x £(i) + D x . (13) 

Proof. This follows from explicit computation. In fact, 

D x = d x + u a k+1 d k a ; (14) 
recalling (11) we have immediately 

[Y {i) ,D x ] = (v? j)ifc+1 - A^), fe ) d h a - (D x Z (i) )d x . (15) 

using (12), the r.h.s. of this can be rewritten as 

-Ufc+lArfo + <Tij(ll>V),k - U k+1 £(*))] 9 a - (D x £ {i) )d x 

= ~{D x $,(i))D x + cry (y (i ) - £(j)d x ) - o-ij (D x - d x ) 
= - + cry D x + CTyYy) ; 

this completes the proof. A 

Remark 2.3 It follows easily from (15) that, conversely, if the Y^ satisfy (13), then 
they satisfy (12), i.e. are a cr-prolonged set. 

Theorem 1. Let y be a set of a -prolonged vector fields, and let be independent 
common differential invariants of order k for all of them. Then 

e - (16 » 

is a common differential invariant of order k + 1 for all of them. 
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Proof. It is obvious that 6 is of order k + 1. To show that it is invariant under any 
of the Y(j), we just proceed by straightforward computation; first of all, by 

[Y (i) (D x Q] • (D xV ) - (AcQ • [Y(i)(D x v)] X 



it is clear we just have to show that the numerator \ vanishes. On the other hand, we 
have (recalling that by assumption Yu\(C) = Yu\{r}) = 0, and using Lemma 1 above) 

X = [Y (l) (D x C)} ■ (D x rj) - {D X Q ■ [Y {i) (D xV )} 

= ([D x (Y {t) C)} + \Y {i) ,D x ]C)-{D x ri) 

- (D x () ■ ([D x (Y (l)V )} + [Y {i) ,D x ] V ) 

= {[Y {i) ,D x }() ■ (D xV ) - (D X Q ■ ([Y {i) ,D x } V ) 

= - (D X Q (D xV ) 

+ (D x () (DzSw+ayfo) (D xV ) = . 

This shows indeed x = and hence the Theorem. A 

Corollary 1. If a complete basis of (independent) invariants of order zero and one 
for y are known, we can successively generate a basis for invariants of all orders. 

Proof. One should only check that the independence of the rj a and (p (differential 
invariants of order zero and k) implies independence of derived invariants. This only 
concerns functional independence of the Qj g := (D x £p/ D x rj a ), and is shown e.g. in 
Olver [23] . Note here (and there) we are using the fact there is only one independent 
variable. A 



Note that in Theorem 1 the aij are arbitrary (smooth) functions. It is clear that 
the relations between the Xu\ will not be shared by the Yu\ for arbitrary cry; the 
condition under which the Y^ have the same commutation properties as the Xu\ (i.e. 
the commutator of cr-prolonged vector fields in y is the same as the cr-prolongation of 
the commutator of vector fields in X), is discussed in Section 3 (under the simplifying 
assumptions that the vector fields do not act on the independent variable). 



3. Involutivity of cr-prolonged sets of vector fields 

The IBDP for cr-prolonged vector fields, and hence the reduction procedure to be 
described in the following, is based on the condition that the cr-prolonged vector fields 
are in involution. In this section we discuss conditions guaranteeing this is the case, 
based on the preliminary condition that the {Xt{\} are in involution. 

Theorem 2. Assume the vector fields = t;(i)d x + f%\d a (i = l,...,n) are in 

involution, with [X^,X^] — jj,^ X^) f or iHj smooth functions on M. Then their 
a -prolongations Y^ satisfy the same involution relations, i.e. \Y^,Y^] = fjb^Y^k), if 
and only if the a/ satisfy, for all k = 1, n, the equations 
{[Y(a 1 k )-Y 1 (a l k )] + 
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Proof. This follows from an explicit computation, described hereafter; we proceed 
by induction on the order of the prolongation. Denoting by Zi the (q — l)-th a- 
prolongation of the X i: we have 

[Y {i) ,Y {j) ] = [Z (i) ,Z (j) ] + [Y (i) (r (jU ) ~ %(Vfi), 9 )] dl := [Z (l) ,Z (j) ] + F^ q d* . 

Thus, assuming [Z^,Z^] = n^Z^ (i.e. the involution relations are satisfied for 
(q — l)-th prolongations), the requirement that [Y^, Y^} — n k jY k is equivalent to the 
requirement that 

*ii, q = 4 r {kU • (is) 

The Fj" needs to be rewritten for easier comparison with the r.h.s. of (18). In 
the course of this computation we will write, for a short notation, for tp^ q _ x ; note 
this entails Y^) = Z,(^); with this, (18) reads 

Ff jiq = ii% (D x n + a km $ a J . (19) 

Using also Lemma 1, with standard algebra we get (omitting the vector indices a for 
ease of writing - and reading) 

= Yi[(D x ^j - u k D x £j) + a jk ($ k - WfeCfe)] 

- Yj[(D x $i - u k D x ii) + o lk {<S> k - u k £ k )] 

= (Y i D x $ j -Y j D x <f> i ) + [Ki(o-jfc)-ij-(^ifc)](*fc-«fc&) 

- [KkYji^k) ~ X jk Yi($ k )] 

= D x [Y i ($ j )-Y j (<S> i )} + [<JikY k (*j) - <T jk Y k ($i)] 

+ [Yi(a jk ) -Yj(cr ik )\$ k - [<J ik Yj{^ k ) -a jk Yi($ k )] 
= D x [Zi{9j) - Zj^i)] + [a ik Z k {<S>j) - a jk Z k ($i)] 

+ [Yt((T 3 k) - Yj{a ik )] $ fe - [o- ik Zj{$ k ) - o- jk Zi(<f> k )} 
= D x ((i%<5> k ) + [Yi(a jk ) -Yjio-ik^^k 

+ <j ik [Z k {^j)-Zj{^ k )] - o- jk [Z k (<f>i) - Zi(<$> k )} 
= (£> x /i£-)*fc + fiijDx^k) + [Yi(<r jk ) -lj(o- ifc )]*fc 
+ <7ikf4j®e ~ °~3kAi$i ■ 
Comparing this with (19), we must require 

F i:j := (D x ntj)$ k + rfjD x ($ k ) + [Yi(a jk ) -YJfofc)]** 

+ a ik n k j<S>i - <T jk fi e ki <f> e = fi%(D x ^> k ) + [i%<j km <$> m . 
That is, eliminating equal terms on both sides and renaming the summation indices, 
(D x ^ j )^ k + [Yi{o- jk )-Yj{o- ik )]<S> k 

+ (J im ll k m fik - CFjm^ni^k ~ V?jVmk$k = 0. 

We can now collect the $ k term, and finally get (reinserting the index a) 

{[Y(o- 3k )-Y 3 (a lk )} + (D x tfj) + a im ^ k mj ~ Ojm vLi - Vij °m k } $1 = . 

We should now recall that := tp? k ^ q _ 1 , and that (20) was the condition under 
which the commutation properties holding up to order q — 1 are also holding at order 
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On the other hand, since by Definition 1 the are functions on J 1 M, and 
D x /J,ij are also functions on J 1 M since are functions on M, all terms within the 
curly brackets in (20) only depends on J 1 M. That is, if (20) is satisfied for first 
prolongations (i.e. for q = 1), it will be automatically satisfied for all higher order 
prolongations (i.e. for q > 1) as well. 

For q — 1 we have = tp^ q _ 1 — ip^.y hence (20) reduces precisely to (17), and 
the proof is now complete. A 

Corollary 2. //, in the hypotheses and with the notation of Theorem 2, it is 

Y,{o k ) = [(D x 4) + o-rvij ~ <x,"V™ " V%*m] . (21) 
then the {Y^} satisfy the same involution relations as the {X^}. 

Proof. In this case equation (17) is automatically satisfied. We stress (21) is a 
sufficient but not necessary condition for (17) to hold. A 

Remark 3.1. The equation (21) actually involves only the first cr-prolongations, the 
cry being functions on J X M. Thus the r.h.s. of the relation (21) is rewritten as 



da, k . ,„ „ t a s da, 



k 



- ^ft}T + ( D ^h+°Me))-d^- ( 22 ) 

Note that the (21) are nonlinear in the a; thus determining suitable a preserving the 
involution properties of a given set of vector fields (under cr-prolongation) is in general 
a nontrivial task. 

Remark 3.2. If the commute, so that ji k j — 0, then (17) reduces simply to 
[Y^-Y^)} Qfa = 0, and (21) to 

Ua^-Yjiaf) = 0. (23) 



Remark 3.3. We stress that the preservation of involution relations is not required 
by the definition of cr-prolonged sets of vector fields, nor it will be required for the 
reduction of (systems of) differential equations based on cr-symmetries. On the other 
hand, the Yu^ should be an involution system, and this is not guaranteed apriori by 
the fact the X^ are in involution, unless equation (21), or at least (17), is satisfied. 




4. Reduction and er-symmetries 

We will now consider the case where a given differential equation, or rather system 
of differential equations, is invariant under a set of a-prolonged vector fields 3^ = 
{Y(i), Y~(jv)}. We will denote, as above, by X = {-X'(i), ■■■,X( N - ) } the set of vector 
fields in M which, upon cr-prolongation, yield the set y. 

We will start by considering a (rather special) situation where the construction 
needed for reduction is particularly transparent, and then proceed to consider 
generalizations; we trust this way of proceeding will help the reader. 
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4-1. Full reduction 

Let us consider a system of ODEs £ = {E 1 = 0, E m = 0} of order q for dependent 
variables u — (u 1 , u n ), i.e. written as F h (x, u, u^) — 0, where h = 1, m. 

We recall preliminarily that if £ is invariant under a set T- of vector fields in 
J q M, then it is equivalent to an equation which can be written in terms of the joint 
differential invariants (of order up to q) of y, see e.g. [4]; we will thus deal with 
equations written in this form. (One also says that the equation admitting y as a 
system of symmetries are equivalent to equations admitting it as a system of strong 
symmetries. We will systematically use the equivalent form, also in Examples below.) 

Theorem 3. Let X be an involution system of rank n of vector fields over M = 
R x R"; and let y be set, also of rank n, of their a -prolongations. Let the system of 
m ordinary differential equations £(x,u^) of order q > 1 in n dependent variables be 
invariant under the set y, i.e. admit X as a -symmetries. Then £ can be reduced to a 
system of m differential equations of order q — 1 . 

Proof. Let us consider the set X of vector fields in M which, upon er-prolongati- 
on, yield the set y. We have supposed the number of dependent variables is the 
same as the rank of the set of vector fields in the set y, and therefore also in the set 
X. This means that the distribution V(X) spanned by the X in TM has dimension 
n, and hence that we can introduce local coordinates (y; w 1 , ...,w n ) in M such that 
Xu)(y) = for all i. (Note this will in general mix x and the u, i.e. dependent and 
independent coordinates; this is analogous to the standard situation [23, 27].) 

Consider now the action of y in J X M, which has dimension d\ = (fin +1); 
again the distribution generated by the y has dimension n, and hence there exist n 
differential invariants of order one (in addition to r\ = y, invariant of order zero). We 
denote these by ...£„). 

When considering the action of y in J q M, we note the latter has dimension 
d q = [(q+l)n+l], while T>(y) still has dimension n, hence we have (qn+1) invariants; 
one of these is rj, and we have n differential invariants of each order. 

The latter can be built using Theorem 1; in the (y;w) coordinates the higher 
order differential invariants (i^ of order k+ 1 < q (thus < k < q — 1) will simply be 

(i k) = D k y ( a , (o=l,...,n) . (24) 

As the system is supposed to be invariant under the y, the differential equations 
can be written in terms of the common differential invariants of y, i.e. in view of our 
discussion as 

F h (y;C,<: il) ,-,C (q - 1) ) = • (25) 

Introducing now new coordinates Zi = Q (i = 1, n), and recalling (24) - which 
implies that (i k ^ = Dy( a = DyZ a - the system (25) will be written as 

F h (y;z,zW,...,z^) = , (26) 

where of course z^ denotes the set of fc-th order derivatives (d k Zi/dy k ). 

We have thus constructively obtained the reduction of the initial set of equations 
F h (x, u, u^) = to a set £ of equations (26) of order (q — 1). 

Finally we note that the requirement q > 1 guarantees that all invariants of order 
q are built recursively from those of order zero and one (see Corollary 1); we could 
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not say the same for q = 1 (as also obvious from the fact the action at order zero has 
no trace of the functions). A 

Remark 4.1. Needless to say, the situation we are considering in Theorem 3, i.e. as 
many symmetries (or more precisely the rank of their system) as dependent variables 
in the system, is highly special and rarely met in practice. As anticipated at the 
beginning of the section, we are discussing this case first as it helps in visualizing 
the mechanism under our general reduction results, which will hold under much less 
exacting conditions. The geometrical meaning of Theorem 3 is quite transparent and 
shows easily the direction of generalization: indeed we are using the fact the system y 
on J q M admits an invariant distribution (of dimension n), and hence common integral 
manifolds, in T(J q M), and use this information to reduce the equations. The fact n 
is precisely the number of dependent variables allows to implement this reduction 
by lowering the order of the full system, while in general we will have less uniform 
reduction, as discussed below. 

Remark 4.2. Note that if we have solutions of the reduced equations £ , i.e. explicit 
functions Zi(y), in order to have solutions of the initial equations £ we must solve the 
reconstruction equations system 

Zi(y) = Ct(y,w,w y ) (27) 

as a differential equations system for the w(y) — {w\(y), w n (y)}. After this, we 
will also have to invert the change of coordinates (x, u) — > (y, w) in order to express 
the solution in the original coordinates in which the £ where set. 

Remark 4.3. We stress that the reconstruction equations (27) are a (generally) 
nonlinear system of non-autonomous equations. Thus, we are definitely not guaranteed 
to be able to solve them. 

Remark 4.4. As the differential invariants thus built are common to all the vector 
fields, they are in particular differential invariants for any one of them. Suppose one of 
them, say Y( ) represents the differential equation we want to reduce. The differential 
invariants can then be used to perform the reduction along the lines of standard 
reduction for ODEs with symmetry, see again Olver [23]. Note also that if we wish 
V(o) to be the standard prolongation of A( ), we should just require that agj = for 
all j. 



4-2. Partial reduction 

In Theorem 3, the set y was of the same rank n as the number of dependent variables. 
Needless to say, this is a very special and fortunate case; in general, we will have (or 
be able to determine) a set of rank r < n. 

It should be clear by inspection that the approach followed in the proof of Theorem 
3 can be followed also in this case, leading of course to somewhat different conclusions. 
We state this as a variant to Theorem 3: 

Theorem 4. Let the system of m ordinary differential equations £(x,u^) in n 
dependent variables be invariant under the set y, of rank r < n, of a -prolonged vector 
fields. Then £ can be reduced to a system of m differential equations, depending on 
derivatives of order up to q > 1 for (n — r) variables and on derivatives of order up 
to (q — 1) for r variables. In particular, if the Jacobian (SF/Su^) is nonsingular, it 
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can be reduced to a system of m equations, r of them of order q — 1, and (m — r) of 
them still of order q. 

Proof. Proceeding as in the proof to Theorem 3, we determine r independent invariant 
functions Q(x, u, u x ), i — 1, ...,r; up to relabelling of variables we can then pass to 
new variables (vi,Zj) defined by 

Vi = Ui for 1 < i < s = (n — r) , Zj = Cj f° r 1 < i < r . 

The functions F h defining the ODEs under study can then be written as 

In other words the system is of order q in the variables v and of order (q — 1) in the 
variables Zf, we can thus rewrite it as in the statement. A 

Remark 4.5. Needless to say, one could as well relax the assumption that all 
equations are of the same order, obtaining a corresponding special result by the same 
procedure. The necessary modifications are clear and we leave details to readers, 
possibly with a specific application in view. 



5. Geometrical features of er-prolongations 

The result of Theorem 1 can appear rather surprising: in fact, the "invariants by 
differentiation" property is intimately related to the standard prolongation structure, 
and one should think that if it works for a somewhat arbitrary (or apparently so at 
first sight) modified structure of vector fields there is some geometrical reason. This 
is indeed the case, as we discuss in the present section. 

Let us consider a set of cr-prolonged vector fields {Yi} in J n M. As remarked 
before (see Section 1) following Pucci and Saccomandi, the invariants only depend on 
the distribution spanned by these vector fields; in other words, passing to a different 
set of vector fields {Zi} which are point- wise in the linear span of the Yi, we have the 
same distribution and hence in particular the same set of invariant functions in J n M. 
That is invariants are related to the module (over C°°(J n M, R) function) generated 
by {Yi}, and are the same for any choice of the module generators. 

We want to discuss if the same module with ^-prolonged generators {Yi} also 
admits generators {Zi} which are standardly prolonged vector fields. We will look for 
these in the form 

Zi = Aij Yj , (28) 

where the are smooth functions on M, and A is non-singular at all points; note 
that as A is a function on M, the restriction of Zi to J h M is a proper vector field on 
J h M (since this is the case for Yi). It clearly suffices to discuss the situation for first 
prolongations, as fc-th prolongations are obtained as first prolongations of (k — l)-th 
prolongations. 

Theorem 5. Let {Xi} be vector fields in M, and {Yi} be their a -prolongations. 
The module generated by the set of a-prolonged vector fields Yi obtained for a = 
coincides with the module generated by the standard prolongations Zi of the vector 
fields Wi = AijXj with A^ smooth functions on M satisfying 



D X A = A a 



(29) 
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Proof. The argument is specially transparent if we deal with vertical vector fields 
in M, i.e. with vector fields acting only on dependent variables, and we will start by 
considering this case. We will work in local coordinates, and write 

* = * = + + ^ ■ (30) 

We look for vector fields Wi of the form 

Wi = Aij Xi = (A v ^) A ; (31) 
the standard prolongation of these will be 

Z * = ^ + [P*^M + A tJ (D x ^)] A . (32) 

In view of (30), (31) and (32), it is clear that (28) is satisfied if and only if 

(D x Aij) ipj + A t3 {D x rf) = (D x <p* + o-jwl) ; (33) 

this is readily simplified getting rid of the AD x tp factors, yielding (after rearrangement 
of summation indices) 

{D x A lk ) <p a k = A i:j a jk <p% ■ (34) 

We therefore get that this is true for any choice of the ip% provided the matrix function 
A satisfies (29). 

The computation is just slightly more involved in the case of general vector fields. 
In this case (30) is replaced by 

Xi = &d x + tfd a , Yi = X t + tfdZ (35) 

with 

C = (D x ^-u a x D x ^) + a^-u^); 

and 

Wi = Xi d x + vtda , Zi = Wi + O?0 x a 
where (standard prolongation formula) 

6? - D xV ? - u a x D xXt . 

Setting Wi — AijXj yields 

Xi = Aijij , rg=A 4j <f% . 

Now, requiring that Zj = AijYj amounts to requiring that 0" = Aijip"; in view of the 
above formulas, this is written as 

AijiD^-ulD^j) + (D x A ik )(ip a k -u a J k ) = Aij[(D x ^-u a x ^) + aj k (<p a k -u a x ^)} ; 

eliminating the equal terms on both sides, this yields 

{D x A ik - AijOj k ) ((ft - <£ fe ) = . 

Again, this holds for any (<p k ,£k) provided A satisfies (29). A 
Remark 5.1. As A is invertible we can also write the relation (29) in the form 

a = A' 1 D X A . (36) 
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Note that the solution to equation (29), equivalently to (36), will in general not be 
unique; see Example 4. 

Remark 5.2. It should be stressed that, as also appearing from (36), the function A 
can be a nonlocal one; in particular, 



exp 



dx 



(37) 



and unless a = D x S(x, u) for some local function S we get a nonlocal function. 

Remark 5.3. It should be stressed that the set {Zi} and the set {Yi} will in 
general not have the same invariants. Moreover, they could have different involution 
properties, see Example 5. 

Theorem 5 has a converse, which we state as a Lemma (Lemma 2) in view of its 
interest to build concrete examples. We also stress that this construction, if applied 
to vector fields which are in involution, will produce vector fields which are again in 
involution (Lemma 3). 

Lemma 2. Let {Wi} be a set of vector fields on M, and {Zi} their standard 
prolongations. Consider linear combinations of these with a nonsingular matrix 
function A : M — > Mat(n), Yi = AijZj. Then the {Yi} are a set of a -prolonged 
vector fields with a -prolongation coefficients o~ij = A (D X A). 

Proof. Obvious by construction. A 

Lemma 3. Let {Yi} be a set of vector fields on J k M, and assume they are in 
involution, so that [Y i} Yj] = H k jY k for some functions : J k M — > R. Consider 
linear combinations of these with a nonsingular matrix function A : M — > Mat(n), 
Zi = AijYj . Then the {Zi} are in involution, [Zi, Zj] = d^Z^, with 8^ : J k M — > R 
explicitly given by = [A iTO /^A^ + A im Y m (A jh ) - A jm Y m (A ih )] AjJ:. 

Proof. For Zi as in the statement, it follows by standard algebra that 

[Zi, Zj] = Ai m Ajk [Y m , Yfe] + [(Ai m Y m (Ajk)\Yk — [AjkY m (Ai m )]Y m 
= [Aim Aj k [i m i~]Yi + [(Ai m Y m (Ajk )]Yk — [AjkY m (Ai m )]Y m 

= [(^imMm^ji + Ai m Y m (^Aj\A — AjmYjxJyAih)} ^hk \ ■ 

This proves the Zi are in involution, and moreover provides the explicit form for 
the functions d\ '■ in terms of the functions u k 4 and of the transformation matrix A, 
invertible by assumption. A 

Remark 5.4. Theorem 5 and its converse (Lemma 2) can also be stated as the 
commutativity of the diagram prol 

A 

x a ^ 01 y 

(where of course X = {Xi}, y = {Yi}, and so on) if A and a arc related by 
g - A (D X A- 1 ) = - (D X A) A^ 1 . 
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Remark 5.5. Suppose we operate with a second (invertible) linear transformation 
B, so that the vector fields Xi — f^d a and Yi = Xi + ip^d^ are mapped into 
Xi = BjXj = ffda and Y = B/Yj = X t + fyd x a . Then 

i>i = B^j = B^D^j+a^ifik) ; 

writing ip — B^ 1 ^, we obtain readily that the vector fields {Yi} are the tr-prolongation 
of the {Xi} (that is, satisfy tpi = D x (pi + a^ipj) with 

a = BaB- 1 + B (D X B~ X ) . 

Note this reproduces (36) if we start from a — 0. Needless to say, this states that cr 
transform as gauge coefficients [11]. 

Remark 5.6. At first sight one could think that Theorem 5 and Lemma 2 state the 
triviality of cr-prolongations, in the sense that sets of cr-prolonged vector fields can be 
mapped into (and hence are equivalent to) sets of standardly prolonged ones. However, 
a little reflection shows that this equivalence is in general only local; in particular if 
the distribution spanned by the Yi is singular at some points (we recall this can mean 
either that some vector fields become singular, or that the rank becomes smaller, at 
specific points) so that the domain on which our procedures are well defined is not 
simply connected, then some cohomological effects will appear, and the equivalence 
will hold only on contractible subsets of M and J h M. The situation here is quite 
similar to the one met in discussing ^-prolongations [8] (for other similarities - and 
differences - between a and //-prolongations, see the Appendix A). 

Remark 5.7. It should be stressed again that the results of this section are a 
direct generalization of those obtained by Pucci and Saccomandi [25] for standard 
A-prolongations. In fact, they showed that - in the scalar ODE case - the most 
general class of vector fields in J k M having the same characteristics as a standardly 
prolonged one (and as a consequence, sharing with them the IBDP) is precisely that 
of A-prolonged vector fields (except for a degenerate case, p\ = p2 = in their 
classification, corresponding to contact symmetries) . Here we have shown that - in the 
case of vector ODEs - the most general class of sets of vector fields in J k M having the 
same integral manifolds as a set of standardly prolonged ones (and as a consequence, 
sharing with them the IBDP) is precisely that of cr-prolonged sets of vector fields. 
If we agree that our technique provides the proper generalization of their "telescopic 
vector fields" (which were defined as those having the same integral lines as some 
standardly prolonged vector field) to the case where one has a system vector fields in 
involution (rather than single ones) and correspondingly looks at integral manifolds 
(rather than integral lines), our discussion confirms their statement that "telescopic 
vector fields seem to be the natural framework for the study of reduction methods 
based on differential invariants" (see [25], p. 6151). 

6. Discussion and possible generalizations 

We have introduced a new modified prolongation operator, called cr-prolongation or 
joint A-prolongation; this does not apply to vector fields individually, but instead to 
a given set of vector fields in involution. We have shown that the cr-prolonged vector 
fields still possess the "invariant by differentiation property" , and hence if they have 
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suitable relations with the vector field describing a system of ODEs they can be used 
for reducing the ODE system. 

We have also discussed the geometrical meaning of er-prolonged vector fields; we 
found that these generalize to higher dimensions the property of A-prolongations - 
discovered by Pucci and Saccomandi [25] - of having the same integral lines as some 
other standardly prolonged vector field, and being the most general vector fields in 
J k M (projectable to J h M for all < h < m) with this property; needless to say, as 
we are in higher dimension integral lines are here replaced by integral manifolds. 

Sets of vector fields which, after being cr-prolongcd (with a specific matrix a), 
leave invariant a given system of ODEs have been christened cr-symmetries. In general, 
vector fields which are cr-symmetries of a system are not ordinary symmetries as well; 
thus our construction really gives new weapons to the arsenal of useful and structurally 
interesting procedures for studying nonlinear systems. 

As already mentioned, the same procedure is studied from a slightly different 
perspective and with some difference in a companion paper [9]. In particular, in that 
paper one is allowing vector fields which are not symmetries of the evolution equation 
but which generate an involution system with the vector field representing it; and 
one also allows vector fields in J k M which are not necessarily prolongations of vector 
fields in M, in the same way as the functions cry are not necessarily depending only 
on variables in J M. From the present point of view, these represent generalized 
a-symmetries (in the same sense as one usually speaks of generalized symmetries in 
standard symmetry analysis [1, 6, 13, 23, 24, 27]); discussing these would be outside 
the limits of the present work. 

The careful reader has probably noted that throughout the paper we discussed 
how to use cr-symmetries if we know them for a given equation, but avoided to discuss 
how to determine these. The reason for this is that the determining equations are in 
this case (as is also the case for A-symmetries) functional equations - and in this case 
matrix functional equations; fully solving them is thus in general far beyond reach, and 
one has to rely on educated guesswork (or sheer luck) in order to determine convenient 
special solutions (see also Appendix B in this respect), i.e. assume a given functional 
form for the cr. Luckily, one does not need to know the most general a-symmetries 
to be able to use them, and special solutions can be enough to reduce the equations 
under study (this feature is again in common with A-symmetries). 

Note that on the other hand, determining the equations or systems which admit 
a given set of vector fields as a cr-symmetry (with assigned a) is a more tractable 
problem; actually, once we have determined differential invariants of order zero and 
one, the IBDP makes it immediate to determine differential invariants of higher orders 
and hence invariants equations and systems. 

Some other directions of further research can be foreseen, and we very sketchily 
discuss them before passing to discuss a number of Examples illustrating our results; 
these are left for future investigations, by ourselves or by some readers of this note. 

(1) Here we mainly discussed the general case, i.e. equations of arbitrary order q > 1 
with n dependent variables and an involution system of rank d < n. The case of 
dynamical systems, i.e. q = 1, would be of obvious interest, and deserves further 
investigation along the lines of the present paper; as already mentioned, it is 
considered from a slightly different point of view in the companion paper [9] . A 
forthcoming paper [10] contains a study along the lines of the present approach. 

(2) Due to obvious physical reasons, one is specially interested in systems of ODEs 
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which arise from a variational principle - i.e. systems of Euler-Lagrange 
equations. These will present special features, and it is an easy guess that a- 
symmetries will be specially effective in symmetry reducing this kind of systems, 
exactly as it happens for standard symmetries (Noether theorem) and for A- and 
^-symmetries [7, 21]. 

(3) Lambda-symmetries are naturally related to nonlocal (standard) symmetries 
[5, 22, 18]. It is natural to expect that, as cr-symmetrics are a generalization 
of A-symmetries, some relations exists between tr-symmetries and some type of 
nonlocal symmetries (possibly based on the construction of Section 5). 

(4) Finally, we would like to mention that an extension of the classical symmetry 
reduction is based on so called solvable structures, see e.g. [2, 3, 12, 14, 26] for 
details. We expect a corresponding generalization, along the lines discussed here, 
would be possible; and actually quite natural as that theory is naturally set in 
terms of distributions of vector fields rather than of single ones. In this respect 
it should be noted that already Pucci and Saccomandi suggested that "a more 
geometric theory of telescopic vector fields and A-symmetries is surely possible by 
means of the theory of solvable structures or the theory of coverings" (see [25], 
p. 6154), so we are again suggesting a generalization of their approach. 

7. Examples 

Example 1. 

Let us consider X = R with coordinate x, U = R 2 with coordinates (u, v). In 
M = X x U the system 



J u xx = u x v x (l + e ") 
\ Vxx = u x v x (l - e~ v ) 
is cr-symmetric under the vector fields 

Xi = d u , X 2 = d v 

with the functions Oij given in matrix form by 



Note these are in involution, and actually commute, [Yi, Y 2 ] = 0. Note also, in passing, 
that in this case equation (21) is satisfied, as easily checked; thus we are in the situation 
claimed by Theorem 2. 

The only common geometrical invariant (differential invariant of order zero) for 
these vector fields is obviously r\ = t; common differential invariants of order one for 
y = {Y\, Y 2 } are provided by 




The cr-prolonged vector fields Y t in J 2 M are then given by 

Yi = d u + v x d Vx + u x v x <9 Uxx + v xx <% xx ; 
Y 2 = d v + u x d Ux + u xx d Uxx + u x v x <9„ xx . 



d 0) = c-« v x , d 0) - e- u x 
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According to our Theorem 1, 
order two. In fact, we have 



(i) _ 



= -D K [Ci ^ should be differential invariants of 



D x 



AO) 



e i^xx ^x v x) 



it is immediate to check that these are indeed invariant under both Y\ and Y 2 , as 
claimed by Theorem 1. 

Let us now consider, to give an illustration of Theorem 3, any second order 
differential equation(s) invariant under 3^; these are necessarily written in the form 

F h (t, zi,z 2 , wi, w 2 ) = , 

where we have set, for ease of writing, 

zi = d 0) , h = d 0) ; = d 1} , - 2 = <$> . 

The above system corresponds to the choice 

F 1 — w\ + ziz 2 , F 2 = w 2 - ziz 2 . 
Through the change of coordinates 

z\ = e~ u v x , z 2 = e~ v u x 

system (38) will reduce, eliminating the common exponential factors, to the first order 
system 

dz\ /dx = — z\ z 2 , 
dz 2 /dx = z\ z 2 . 

Example 2. 



(39) 



Let us again consider X = R with coordinate x and U — R 2 with coordinates (u, v). 
We take now the vector fields 



X x = ud u , X 2 



-ud v ; 



these are in involution, and actually satisfy [Xl,^] = X 2 . 
We consider the functions given in matrix form by 



a = 



u 
u T 



The a- prolonged vector fields Yi in J 2 M are then given by 
Yi = u d u + u x d Ux + u 2 d Vx 

-\- (ll U x H~ l^xx ) ^«xx + (3wwx) <9„ xx ; 
Y 2 = - ud v - uu x d Ux + u x d Vm 

- (uu xx + 2u 2 x ) <9 Uxx - (u xx + u 2 u x ) d Vx 



d 1} = d x 



AO) 

^2 
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Note these are in involution, and actually satisfy [Yi,Y" 2 ] = Y 2 ; these are the same 
involution relations satisfied by Xi and X 2 , and also in this case one can check that 
equation (21) is indeed satisfied (and thus Theorem 2 holds). 

The only common geometrical invariant (differential invariant of order zero) for 
X = {X\,X 2 } is obviously r\ — x; common differential invariants of order one for 
y = {Y\,Y 2 } are provided by 

d 0) - e""^, d 0) = 2^- u 2 -2(l-e-)^. 

According to our Theorem 1, — D x [(^] should be differential invariants of 
order two. In fact, we have 

C[ 1} = D x Ci° } = (e~ v /u 2 ) (uu xx -u 2 x - uu x v x ) , 

(2/x 2 ) [{u 2 x - u 3 u x - uu xx + u 2 v xx ) 

+ (uu xx -u 2 x - uu x v x )] . 

It is immediate to check that these are indeed invariant under both Y\ and Y 2 . 
Let us now consider, with a view at Theorem 3, the second order equations 

u xx = (l/u) [(4e~ v - 7)u 2 x + (9uv x - 4u 3 )u x 

+ e v (m 6 + 4(m 2 + u 2 v 2 x + v?u x - u^v x - 2uu x v x ))] , 

v xx = {l/u 2 ) [u 3 u x + (e~ 2v - l)u 2 x + uu xx 

- e~ v (uu xx - 2uu x v x + (l/2)u 3 u x )] . (40) 

These are invariant under both Y\ and Y 2 , as can be checked by explicit 
computations; passing to the F-adapted coordinates Ci, C2 requires to determine the 
inverse change of coordinates, which is (writing Wi = dQ/dx) 

u x = «Ci , u xx = (w/2)[2e t 'u;i+Ci(u 2 + 2Ci + C2)] , 

v x - (u 2 + ( 2 )/2 , v xx = {l/2)[2(e v - l)wi + w 2 + (3u 2 + C 2 )Ci] • 

Plugging these into the equations (40), these reduce to 

wi:=d(i/dx = C|, w 2 :=d( 2 /dx = Ci C2 , (41) 
i.e. to a system of first order equations. 

Example 3. 

In the previous two examples we have just remarked equation (21) was satisfied; in 
order to show this is not always the case, and that failing to satisfy this equations 
leads in general to er-prolonged vector fields which are not in involution, we consider 
again the vector fields of Examples 1 and 2 but with different matrices a. 

Case 1. In the case of Example 1, i.e. X\ — d u and X 2 — d v , and a matrix a which 
is just the transpose of that used in there, i.e. 



a = 



u x 
v x 
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In this case the first joint prolongation is given by Y 1 = d u + u x d Vx , Y 2 = d v + v x d U:c ; 
thus the involution relations between the X arc not mapped into corresponding 
relations between the Y; more precisely, [Yi, Y 2 ] = Y 3 := u x d Ux — v x d Vx . Moreover, in 
order to close the involution relations we also have to introduce two other new vector 
fields, Y4 u x d Vx and Y5 := v x d Ux . With these, the involution (actually, algebraic - 
due to the constant coefficients in X\ and X 2 ) relations are then 

[Yi,Y 2 ] = Y 3 , [Fr, Y 3 ] = -2 y 4 , [Yi, Y 4 ] = , [Y U Y 5 ] = Y 3 , 
[Y 2 , Y 3 ] = 2 Y 5 , [Y 2 , Yi] = -Y 3 , [Y 2 ,Y 5 ] = , 
[Y 3 , Y A ] = 2 F 4 , [Y 3 , Y 5 ] =-2Y 6 , [Y 4 , Y 5 ] = Y 3 . 

Thus, the involution system y is actually an algebra; and the vector field we have 
to add are an ideal in y. 

As for equations (21) and (17), it is immediate to check these are not satisfied: 
the r.h.s. of both vanishes (in this case /i^ = for all i,j,k), and the l.h.s. of (21) 
is easily checked to be nonzero. More precisely, let us define = Yi(aj) — Yj(af) 
(it is obvious that — —Qji, and hence Q\ = for all i); then Q 12 = {u Xl —v x ). 
Applying this on vectors (pk, i-e. computing the l.h.s. of (17), we get again (the matrix 
built with the ipi as columns is just the identity) 

QijVk = («*. - v *) + • 

Case 2. Similar considerations apply if we consider the vector fields of Example 2, 
i.e. Xi = ud u and X 2 = —ud v , and a matrix which is the transpose of the one used 
there, 

° = (° "<;) - ^ 

In this case Y\ — ud u + u x d Ulc — uu x d Vx ; Y 2 = ud v + u 2 d Ux — u x d Vx . We have then to 
introduce a vector field Y 3 = m 3 9 1>x and have 

[Yi,Y 2 ] =Y 2 + Y 3 , [Y U Y 3 ] =3Y 3 , [Y 2 ,Y 3 ] = . 

In this case we get Qi 2 = (u, — u 2 ); Qi 2 f% — (u 2 , u 3 ). 

Case 3. Finally, let us consider the vector fields of Example 1, X\ = d u and 
X 2 = d Vl and the matrix a just considered above, see (42). In this case the first 
joint prolongation is given by Y\ = d u + u x d Vx , Y 2 = d v + ud U!c ; in order to close the 
involution relations we have to add two new vector fields in J 1 M, i.e. Y 3 = d Ux —ud Vx ; 
Y4 = d Vx . With these, the involution (actually, again algebraic) relations are then 

[Y U Y 2 ]=Y 3 , [Y 1 ,Y 3 ] = -2Y 4 , [Fi,F 4 ]=0 , 
[r 2 ,r 3 ] = o , [Y 2 ,Y 4 } = , [Y 3 ,Y 4 }=0 . 

As for equations (21) and (17), again these are not satisfied; the r.h.s. of both 
vanishes for fi^ = 0, and as for the l.h.s. we have Qi 2 = (1, — u), Qijfl = (1, —u) ^ 
(0,0). 

Example 4- 

We aim now at illustrations of Theorem 5. Let us consider the (obviously commuting) 
vector fields X\ = d v , X 2 = (l/u)d u ; we a-prolong them with the (nearly trivial) 
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matrix a — (u x /u) I (I being the two-dimensional identity matrix); in this way we get 
Y 1= d v + {u x /u)d Vx , Y 2 = (l/u)8 u . 

We note that these are not in involution: in fact, [Y^Y^] = ( u x/u' i )d V:c := Y3; 
this satisfies [Yi, Y 3 ] = 0, [Y 2 , Y 3 ] = — (3/u 2 )Y 3 , hence closes the algebra. 

With the notation introduced in Example 3, we have Qi 2 — (u x /u 3 , 0); Q\ 2 ^%. = 
(0,u x /u 3 ). 

We now look for a set of standardly prolonged vector fields which, as stated by 
Theorem 5, are in the module generated by (Y\, Y 2 ). As suggested by Theorem 5, we 
look for A(x 7 u, v) solution to (29), i.e. in this case to 

D X A = {u x /u) A . 

The solution to this equation is not unique: any matrix of the form 



K 



k\u k 2 u 
k 3 u k^u 



with ki constants is a solution. Let us look in particular at the matrices 

ii\ _ (u 

u OJ ' V0 u 

In the present case, using A we get 

Z[ A) = ud u + u x d Ux , Z { 2 A) = d v ; 
using B we get instead 

Z[ B) = ud v + u x d Vm , Z { 2 B) = d u . 
These sets are transformed one into the other by the matrix 



M = 



1 

1 



which also maps A and B one into the other; more precisely, writing 

z\ A) = v ?d a + %&z , z[ B) = 7 ?d a + etdi , 

we have 

We stress that 

[Z[ A \Z ( 2 A) } = 0; [Z[ B \Z { 2 B) ] = -d v # span(^ B) ,^ 



( B) ( B) ( B) 

On the other hand, Z\ = —d v commutes with both Z\ and Z 2 . 

Finally, let us look at invariants; a set of common invariants (beside x) for {Z^} 
is provided by {v x ,u x /u}. As for the vector fields Z\ , we need to add the third 
vector field Z\ ' and hence we expect (as is indeed the case) that only one common 
differential invariant of order one exists (beside x); this is u x . 
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Example 5. 

We now illustrate Lemma 2, i.e. the converse of Theorem 5. Let us consider the vector 
fields W\ = ud u + vd Vl W 2 = —vd u + ud v ; their (standard) prolongations are of course 

Z x = ud u + vd v + u x d Ux + v x d Vx Z 2 = -vd u + ud v - v x d Ux + u x d Vx . 

According to Lemma 2, if we operate on these with a linear transformation 
A{x,u, v) depending on variables in M, we should obtain a set of vector fields which 
are a cr-prolongcd set, with a corresponding to a = A (D X A); more precisely, the 
{Yi} should be the cr-prolongation of the Xi = AijWj. 

We will write $ = (1 — wi;)^ 1 , and choose 

*- (I I) • A ~ l ■ • (T -„ 

which yields 

a = $ (- vu - v * 
\ u x -uv x 

With simple computations, we get 

Yi = $ [(1 - u)vd u + (1 - v)ud v + (v x - vu x )d Ulc + (u x - uv x )d Vx ] , 

Y 2 = $ [(u + v 2 )d u - (u 2 + v)d v + (u x + vv x )d Ux - (v x + uu x )d Vx ] . 

These are indeed the cr-prolongation, with the a given above, of the vector fields 

X 1 =A lj W j = $ [(l-u)vd u + (l-v)ud v ] , 

X 1 =A lj W j - $ [{u + v 2 )d u - (u 2 +v)d v ] . 

Let us now consider first order differential invariants; it is easily checked that Z\ 
admits as invariants {x, u/v, u x /v x }, while for Z 2 we get {x, u 2 + v 2 , u 2 +v 2 }\ common 
differential invariants are 



u 2 x +v 2 x 
u 2 + v 2 



Q = arctan (u x /v x ) — arctan(u/w) . 



Applying Y\ and Y 2 on these, we obtain a non-zero result in all cases, as stated in 
Remark 5.3. 

Example 6. 

We give another illustration of Lemma 2, this time for vector fields which are not 
defined in u = v = 0. We write p = (u 2 + v 2 ) and choose 

Wi = p~ 1 h(x,u)d u , W 2 = p~ 1 h(x,v)d v , 

with h an arbitrary smooth function; this yields as first standard prolongation the 
vector fields 

Z 1 = Wi + p- 2 (p[D x h(x,u)} - h(x,u)(D x p)) d Ux , 
Z 2 = W 2 + p- 2 (p[D x h(x,v)} - h(x,v)(D x p)) d Vx . 
Let us now consider a transformation 

A = — ( v u ) , A-' = - ( v u 
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and the vector fields 

Xi = A/ Wj , Y t = A/ Z j . 

Writing P = h(x,u), Q — h(x,v); and moreover V = (p(D x P) — 2(D x p)P), 
Q = (p(D x Q) — 2(D x p)Q), these are given explicitly by 

X 1 = - p-^vPdu + uQ8 v ] , 

X 2 =p~ 1 [uPd u - vQd v ] ; 

Y x =Xt - p- 2 [vVd Ux - uQd Vx ] , 

Y 2 =X 2 + p- 2 [uVd u=c - vQd Vx ] . 

It is easily checked that the {Yi} are the cr-prolongation of the {Xi} with 

a =- (~( uu * + vv *) uv x -vu x \ = A , DA -i h 
p \ -uv x + vu x -(uu x + vv x )J 

Example 7. 

In the previous example, both the Yi and the Zi fields were singular. Situations where 
the standardly prolonged fields arc singular but the cr-prolonged ones are regular would 
be of interest; here we deal with such a situation, obtained through a small variation 
on the setting of Example 6. We use the same notation introduced there. 
We consider the (singular) vector fields 

X 1 = (P/p) d u , X 2 = (Q/p) d v ; 

these have standard prolongations 

Z 1 = X 1 + [(p(D x P) - P{D xP ))/p 2 ] d Ux , 
Z 2 = X 2 + [(p(D x Q) - Q(D xP ))/p 2 ] d Vx . 

Acting now with 

A = p 2 ( V ~ U ) 

we get the (regular) vector fields 

X x =p [(Pv)d u - [Qu)d v ] , X 2 = p [{Pu)d u + (Qv)d v ] ; 
Y l =X X + v[p{D x P) - P{D x p)]d Ux - u[p(D x Q) - Q{D x p)]d Vx , 
Y 2 =X 2 + u[p{D x P) - P(D x p)]d Ux + v[p(D x Q) - Q{D xP )]d Vx . 
One can check that the {Yi} are the cr-prolongation of the {Xi}, with 

a = A(D X A^ 1 ) = -(D X A)A~ 1 = -- (^{ uu x + vv x ) uv x -vu x \ 

p \-uv x + vu x h{uu x + vv x )J 

Example 8. 

Let us now see an example of the situation considered in Theorem 4; that is, we will 
have a set of three differential equations of second order, admitting as symmetries a 
set of two CT-prolonged vector fields in involution (actually, commuting). (A number of 
trivial examples are also easily obtained by adding to any n-dimcnsional example for 
theorem 3 (with y of rank n) a new equation for a new dependent variable w(x); this 
is not acted upon by, nor entering in, the coefficients of the considered vector fields 
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and hence no reduction is possible on it. We will thus end up with a system of one 
second order and n first order equations.) 
We consider the equations 

E\ := e~^ u+w ^ (u xx - v xx ) (u x - v x - w x ) - [u x + w^ 1 (u - v) x 
x (-u x 2 + u xx +u x v x + v x w x + w x 2 + w xx ) = ; 

E 2 := e ~ { - u ~ v ~ w) (u x - v x ) (u xx + w xx -u 2 x + u x v x + v x w x + w 2 x ) + 
- [u x - v x - w^ 1 (u - v) x 

x (-u x 2 + u xx + u x v x - v xx + v x w x + w 2 - w xx ) = ; 
£ 3 := (u-v) - e-< u +^ (u x - v x ) x 

x (u xx - v xx - w xx - u 2 + u x v x + v x w x + w 2 ) = ; 
and the (autonomous) commuting vector fields 

X\ = d u + d v — d w , X2 — d u + d v . 

As for their (second) cr-prolongation, we set 

(0 u x + w x ~ 

u x - v x - w x 


and hence we get, with standard computations, 

Y\ = (d u +d v - d w ) + (u x + w x ) (d Ux + d Vx ) 

+ [u xx + w xx + (u x -v x - w x )(u x + w x )\ (9 Uxx + a„ xx ) 
- (u x - v x - w x )(u x + w x ) d Wx;i: ; 

Y 2 = {d u + d v ) + (u x -v x - w x ) (d Ux + d Vx - d Wx 

+ [u xx - v xx + (u x -v x - w x )(u x + w x ) - w x x] (<9 Uxx + d Vx J 

{ u xx Vxx V^xx) d Wxx . 

Now, by trivial dimension counting, there must be two invariants of order zero, 
i.e. on M; these are 

77l = X , T]2 — u - v . 

The common differential invariants of order one for these vector turn out to be 

Ci 0) =e- {u - v - w \u x + v x ) , C, { 2 ] = e- {u+w) (u x -v x -w x ) , = u x - v x . 

The second order differential invariants can be readily computed by the IBDP (using 
rji) as (j 1 ^ = D x C,f \ which yields 

Ci 1 ' = e- {u - v ~ w) (u xx + w xx + u x v x + v x w x -u 2 x + w 2 x ) , 
(jf> = e -("+") (u xx - v xx - w xx + u x v x + v x w x -u 2 x + w 2 x ) , 
A 1 ) _ _ 

S3 u xx V xx . 

these are easily checked to be indeed invariant under the Yi. 
We can then perform the change of variables 

t >(0) >(0) 

i = u - v , zi = Ci , z 2 = Q' ; 
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this of course entails a corresponding change for derivatives, which we do not write 
down explicitly. In the new variables, the equations are written as 



Ei 
E 2 
E 3 



£ ((zi)x/zi) - z 2 ixx = ; 
£ ((z 2 )x/z 2 ) - £x (zi) x = ; 
((z 2 )x/z 2 ) - ((si)x/zi) = . 



The system is now of first order in the Zi variables, and of second order in £; it can be 
rewritten as 

^£ = e dzi = _e_ dz2 = ^_ 

dx 2 Q zi z\ ' dx £2 Z2 ' ^ ^ 

Needless to say, the same procedure - at least up to rewriting the system as first 
order in the Zi variables, and second order in £ - would work for any starting system 

of equations in the form Ei := Fi r\\ , r\ 2 ; ([°\ $\ £^ ; = 0, leading to 

equations written in the new variables as Fi[x, £; z\, z 2 ,£ x ; (zi) x , (z 2 ) x , £xx] = 0. 

Example 9. 

Let us now consider, again in order to illustrate Theorem 4, a case with three 
dependent variables (it, v, w) and two (non commuting) independent vector fields in 
involution, 

Xi = ud u + wd w ; X 2 = -ud v ; [X 1 ,X 2 ] = X 2 . 
We will use the (2x2) matrix 

u x 


With this choice, the a-prolongations of the vector fields {Xi,X 2 } to J 2 M are 

Y\ = u d u + wd w + u x d Um - uu x d Vm + w x d Wx 

"i" u xx (2 u\ + UU XX ) <% xx + 

Y 2 = - ud v - u x d Vx - u xx d Vxx . 

These still satisfy the same commutation (hence involution) relation: [Yi,!^] = Y 2 . 

With trivial algebra one finds the common invariants of order up to two for these 
two vector fields; by trivial dimension counting they must be eight, one of them being 
the trivial one Joo = x. 

A simple basis for the (seven) nontrivial ones is given by: Joi = {w/u), Jn — 
(w x /u), J 21 = (wxx/u); J12 = (u x /u), J 22 = (uxx/u); J13 = [v x + (uUx/2) - (u x v/u)], 
J23 = [v X x + u 2 , + (uiixx/2) — (uxxv/u)]. Another possible, maybe slightly less simple 
but more convenient, basis is provided by 

AO) 1 A 1 ) 1 \ / 2 A 1 ) 1 

Ci = w / u ; Ci = \ uw x - u x w)/u , C2 = u x/u , 

(3^ =v x + UUx/2 - u x v/u ; 

C ± 2) = (u 2 w X x - uwuxx - 2uu x w x + 2wul)/u 3 , 

C 2 ] = (uUxx - U 2 x )/u 2 , C3 2) = {Vx + UUx/2 - VUx/u) . 

This is convenient in that (wherever applicable) 

ci i+1) = A4° • (43) 
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Thus any system of equations : [x; u^ 2 \ v^ 2 \ w^] — of order not higher 
than two written in terms of these joint invariants will admit the {Xi,X 2 } as a- 
symmetries (with the a given above) , and conversely. Note all of these (if nontrivial) 
will necessarily be singular for u = 0. 

In order to consider a concrete example, let us look at the system 

2wu 2 + u 2 (2vu x — 2uv x + w x — w) — uu x (u 3 + w + w x ) 



Uxx — 



V, 



uw 
A + B + C 



w 



2u 2 w 
(u + u x ) (w - w x ) 



XX 



U 

.2-, 



A = 2u°v x - ((Av + u )u x - w + w x )u + (u x (w + w x ) - Avv x )vf, 
B = [(-3wul + Av 2 u x + 2v(w x - w)]u 2 , 
C — 2(ww x — u x (v(w + w x ) — v x w))u + 2u 2 vw . 
In terms of the invariants, this system reads simply 

cf - cf ) + d 0) = o , 
cf - d 1} = o , 
d 2) - 2C« = o. 

In order to reduce it to a system of one second order equation and two first order 
ones, we should change (dependent) coordinates. We choose 

C = d 0) , v = C 2 (1) , 9 = <g> ; (44) 

and recall (43). The equations are then readily expressed in the new coordinates, 
providing 

d 2 t _ dp drj dC 

= 2p > T x =71 ^ Tx = Tx~ e; (45) 

i.e., as claimed, a system of two first order and one second order equations. 



Reduction of ODEs via asymmetries 



26 



Appendix A. Relations between a- and /x-prolongations 

The picture emerging from the discussion of Section 5 above, see in particular Theorem 
5 and Lemma 2, is quite reminiscent of the one holding for ^-prolongations [8, 11] (here 
we assume the reader to be familiar with them; see [11] for details and a bibliography 
on twisted prolongations in general): in this appendix we will discuss similarities and 
differences between the two, and their relations when applied to a given set of vector 
fields. 

Al. Similarities and differences 

First of all, let us stress the analogies between a- and ^-prolongations: in both cases, 

(a) the vector fields Fj obtained from X t through the modified (er or /t) prolongation 
operation are the standard prolongations Zi of different vector fields 14^; 

(b) the relation between Xi and Yi on the one hand, and W, and Z, on the other, 
can be summarized by the action of a linear transformation A; 

(c) this A then determines the relevant matrix (a in our case, a single A in the case of 
/i-prolongations for ODEs) entering in the definition of the modified prolongation 
operation - and reducing to zero for standard prolongations. 

An essential difference between the two cases should however be emphasized: in 
the case of /i-prolongations, the linear transformation acts on the basis vectors in TM 
and more generally in T(J k M) [11]; here instead we are changing the generators of 
the module of vector fields - i.e. act on a structure superimposed to T( J k M) with no 
action on the space itself. 

In other words, in the case of /i-prolongations all vector fields on M and on 
J k M are affected by the action of the linear transformation - as obvious since this 
is a point-dependent change of basis, i.e. a gauge transformation [11] - while here 
we are just changing our choice for the generators of a given module of vector fields, 
or equivalently the generators for a given distribution on TM and T(J fe M), with of 
course no action on the general set of vector fields in M or J k M. 

As a consequence, the /i-prolongation operation acts individually on each vector 
field; they are thus well defined, and indeed different from standard ones, also for single 
vector fields. Moreover - in the framework of ODEs, which is the only one of interest 
here - /^-prolongations reduce to A-prolongations only if the linear transformation A 
is a multiple of the identity, the proportionality factor being a smooth function on M. 

On the other hand, the ^-prolongation operation is defined on finite sets of vector 
fields, and reduces to a somewhat trivial one in the case where the set mentioned above 
is made of a single vector field. They again reduce to usual A-prolongations only if the 
linear transformation A is a multiple of the identity, the proportionality factor being 
a smooth function on M. 

We summarize the situation in the following commutative diagram: 

{Xi} A {Wt} A {Xi} 

J^/j.— prol s ["^ r ° 1 J^cr— prol (46) 

{%} A {Zi} A {Y t } 
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Here Yi is the /z-prolongation of X t , the {Yi} are the cr-prolongation of the {X t }, and 
the inertible linear transformations S and M are related to a and to [i — Adx via 

A = M{D x M~ l ), a = SiDxS- 1 ). (47) 

A 2. Relations 

It is quite clear that albeit the two operations are conceptually different in general, a 
relation exists between the two points of view when we consider a given transformation 
on a given set of vector fields; this is quite similar to the relation between "active" 
and "passive" points of view in fluid mechanics. 

We discuss it for vertical vector fields only (they can be thought as evolutionary 
representatives of general vector fields in M, see [1, 6, 13, 23, 24, 27]), and refer to 
(46) for notation. 

Consider vector fields Wi in M and their prolongations Zi in JM; consider also 
a linear (point-dependent) transformation in M, under which Wi are mapped into Xi 
and Zi into Y t ; we know that the {Yi} are then the cr-prolongation of the {Xi}. We 
want to discuss if there is a linear map M such that these particular Wi and Zi are 
mapped into Xi and Yj (each of them being the ^-prolongation of the corresponding 
X^, with the additional properties that (a) Xi = Xi; (b) Yi = Yi. It will suffice 
to work on first prolongations, as prolongations of order k + 1 can be seen as first 
prolongations of prolongations of order k. 

It is convenient to work in coordinates, writing Wi — <pfd a} and embodying the 
components of the different vector fields into a single matrix $ with elements 

= fl • (48) 

Note that if we have r vector fields with n components, the matrix $ is (n x r). In 
particular, it is a square - and thus possibly invertible - matrix if and only if r = n; 
this is precisely the case we have considered in the main body of the present paper. 
Note also that if r = n, the condition for $ to be invertible (at all points) is precisely 
that the vector fields Xi are independent (at all points). 

Coming back to our vector fields, Xi has components <pf = (S/ifj, while Xi has 
components ipf = M a h ^\. Requiring the two to be the same, i.e. Ip = (p, amounts to 
requiring 

Sjtf = M\y\ ; 
or, acting with M _1 from the left, 

^ b = (M -1 ) 6 a <fj = (M-y a ^(s T y t . 

Rearranging the summation indices and using the notation (48), this is also rewritten 
as the matrix equation 

$ = M -1 $ S T . (49) 

It should be recalled that in our discussion of er-prolongations we have supposed 
the vector fields Xi to be as many as the dependent variables, and to be independent 
at all points of M; in view of the requirement that also S and M are invertible at 
all points, this means that the Wi will also be independent, and hence $ will be 
nonsingular, at all points of M. 
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Thus if we are requiring the M which satisfies (49) for given <& and S, the answer 
is that this is given by 

M = $S' T $~ 1 ; (50) 

similarly, if $ and M are given and we ask if there is a S such that (49) is satisfied, 
we get 

S = $ T M T ($" 1 ) T , S T = $ _1 M$. (51) 

Finally, if M and S are given and we wonder which are the vector fields Wi such that 
Xi = Xi, the answer is provided directly by (49). 

Our discussion so far only concerns vector fields in M; let us now consider also 
vector fields in JM. We write Zi as 

Z, = (d/du a ) + c {d/du a x ) ; 

here of course ipf = D x ipf. We will embody the coefficients tjjf into a matrix VP with 
entries = r/>£; this satisfies * = D x $. 

Proceeding in exactly the same way as above, we easily get that 

Y = X t + [S^ (D x ipj)] (d/du x ) , Y = X t + [M a b {D x ip\)\ (d/du a x ) ; 

thus - assuming Xi = Xi - we have to require 

* S T = M * . 

Recalling now that ^ = D x $, this yields with trivial manipulations 

(D X M) = ^(D^)®- 1 ; 

recalling also (50), we arrive at the condition 

[-f- 1 (D x $)} S T = S T [-f- 1 (D x $)\ . (52) 

In other words, the equality Xi = Xi will be lifted to a corresponding equality 
Yi = Yi between prolonged vector fields if and only if (52) is satisfied. 

A3. Combining fi- and a -prolongations 

The previous discussion also suggests how to combine [i- and cr-prolongations into 
a single prolongation operation; we will call it x-P r °l° n 9 a ^ on (where x stands for 
"combined"): working with matrices $ and * = D x $>, and writing R = S T for ease 
of notation, one should pass to 

$ = M~ l ^R, $ = M-^R; (53) 

the operator mapping $ into will be that of the prolongation combining /i and a 
ones. Inverting the first of the relations (53), we have 

$ = M&R- 1 , 

hence the second of (53) reads 

$ = M'^D^M^ R'^R 

= M- 1 [(D X M)^R- 1 + MiD^R- 1 + M^>(D X M- 1 }R 
= (£> x $) + [M" 1 (D X M)] $ - $ [R- 1 (D X R)] 
= (D x %) + A$ - $9 ; 



Reduction of ODEs via asymmetries 



29 



here we have of course defined 

A := M" 1 (D X M) , 6 := JT 1 (D X R) . 

Reintroducing indices and passing to the vector notation, we rewrite the above relation 
as 

ft = (D x tf) + K\ v \ - (6 T ) i >« ; (54) 

this defines the (first, and then recursively higher order) prolongation operator. 

This can be described in terms of commutative diagrams by complementing (not 
substituting!) the diagram (46) with the diagram 

{x,} -A {PJ A to 

J^ju— prol j^x-prol |cr-prol (55) 

{%} -A {Qi} A {Yi} 

One could thus think of investigating x-symmetries of differential equations, 
defined as sets of vector fields which, when x-prolonged, give sets of vector fields 
which leave the solution manifold of a given set of ordinary differential equations 
invariant. 

However, it is well known that while /x-symmetries are "as useful as standard 
symmetries" in the search for special solutions to PDEs or systems of PDEs, and they 
reduce to usual A-symmetries for scalar ODEs - and are thus in a way again as useful 
as standard ones in their reduction - they are of no (known) use in the case of systems 
of ODEs. Thus at the moment we cannot describe any use of such x-prolongations, 
and they remain at this stage a mathematical curiosity. 
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Appendix B. Determining equations for er-symmetries 

In the Conclusions (Sect. 6 above) we have mentioned that the task of determining 
all a-symmetries for a given equation or system is in general well beyond reach. In 
this Appendix wc will illustrate this statement by a concrete example, i.e. the system 
considered in Example 1 above. This reads 

u xx = u x v x (1 + e~") , 

v xx = u x v x (1 + e~ v ) . (56) 

We will look for a set of two vector fields in involution giving a cr-symmetry, and 
in order to keep formulas not too long, we will search for these with = 0. In this 
case the second prolongations will be written as 

Y (i) = tpfod* + r (i) d^ + xU di a X) . 

where we have denoted '■= (d/du x ) and di xx ^ := (d/du xx ); and the coefficients 
entering in first and second prolongations are now denoted as tp and \ to reduce the 
number of indexes. An explicit computation shows that 

xw = + en r U) 

= D 2 xV a (l) + 2a ij (D xV > a {j) ) + [D x< y l0 + pfo . 

In order to obtain the determining equations for er-symmetries we should apply 
both prolonged vector fields to both equations, obtaining in all cases zero upon 
restriction to the solution manifold of the system (56). 

The determining equations - under the present simplifying assumption = 
are thus in this case (here and below i = 1,2, a = 1,2, and S represents the solution 
manifold to the system) 

1.2 



X(i) - [i>li)V x + u x ip 2 {i) ](l + e u ) + {u x v x )ifffae 



0.(57) 



Substituting according to the above for ip and Xi an d making explicit the 
derivative terms, we get (with <p^ b — dip^/du b , etc.; sum over repeated indices 
is implied) for the left hand side of (57) before restriction to S 



+ 



da i3 



dx x du a x 
Using now (56) to restrict on S, we obtain 



+ 2a tj (^ i)iX + u^ i)i6 ) (58) 



do; 



+ u 



(*4u x (l 



• e -«°))^' 



dx x du a v 1 " " du% 

These are the determining equations we were looking for. It should be noted 
that they cannot be solved according to the standard procedure for determining 
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equations of Lie-point symmetries, as they depend on the unknown matrix function 
a = a(x,u,u x ). 

The best one can do is to look for special solutions of these; e.g., if we look for 
solutions such that <p^ are constant, and a only depends on u x and v x , this reduces 



to 



(^(l + e"" ))^ 



+ o- im a. 



. 



(59) 



Despite its innocent-looking shape (due to compact notation) this is still a system of 
four coupled nonlinear PDEs. 

We will now look for a in the form 



a = 





B(u Xl v x ) 



A(u x ,v x ) 




0, 



= 0; 



0, 



0. 



we also write <p^ = Ci, <p^ = c 2 , <p^ 2 ) = fci, <^ 2 ) = k 2 ; in this way the determining 
equations (59) read 

Ae v (-B Cl e u + (1 + e u )(k 2 u x + k lVx )) 
-(A Vx e u (l + e v )h + e v ( Cl + A u=c (1 + e^k^u^ 
Ae u (-Bc 2 e v + (1 + e v )(k 2 u x + k lVx )) 
-(c 2 e u + A Ux (1 + e u )e v k 2 + A v ^(\ + e v )k 2 )u x v x 
Be v {-Ae u k l + (1 + e u ){c 2 u x + Cl v x )) 
-{B VxCl e u {l + e") + e v {B UxCl (l + e u ) + h))u x v x 
Be u (-Ae v k 2 + (1 + e v )(c 2 u x + Cl v x )) 
-{B Ux c 2 {\ + e u )e v + B Vx c 2 e u (l + e v ) + e u k 2 )u x v x 

Determining the general solutions of these is not easy despite the several ansatzes 
considered to simplify them from the original form. A further reduction is obtained e.g. 
looking for A, B as linear functions, A(u Xl v x ) = p\u x +p 2 v Xl B(u x , v x ) = q\u x + q 2 v x . 
In any case, it can be checked that the vector fields X\ , X 2 and the matrix a considered 
in Example 1 satisfies these equations. 

Note that even assuming a = cr(x, u) (i.e. that are functions on M and 
not on the full J X M), albeit in principles one gets a more tractable problem, in 
that all dependencies on derivatives are explicit and the determining equations can 
be decomposed into equation for the vanishing of different monomials in the u x , 
the resulting PDEs arc nonlinear in the cr^, and hence in general cannot be fully 
solved. E.g., in the case of equations (56) even with the simplifying assumptions 
<P(i) = < /'(i)(' lt ! v ) an d = <J ij{ u i v )i after easily determining that must be of the 
form 



X {i) = 



exp[(l + e~ u )v] 
l + e~ u 



P(i){u) + Q(i){u) 



exp[(l + e~ v )u] 
1 + e- v 



R(i)(u) + Si(u) 



d v , 



one remains with PDEs which are nonlinear in the aij(u,v) and appear to be 
untractable. 
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